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to examine possible acausality in this more complex setup. Despite observing significant 
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1. Introduction 

Recently there has been considerable interest in the novel suggestion that we live in a 
Universe that possesses more than four dimensions. The standard model fields are assumed 
to be confined to a hypersurface (or 3-brane) embedded in this higher dimensional space, 
in contrast the gravitational fields propagate through the whole of spacetime [1-10]. In 
order for this to be a phenomenologically relevant model of our universe, standard four- 
dimensional gravity must be recovered on our brane. There are various ways to do this, 
the most obvious being to assume that the extra dimensions transverse to our brane are 



compact. In this case gravity can be recovered on scales larger than the size of the extra 
dimensions [5-7]. This is different from earlier proposals since the restrictions on the size 
of the extra dimensions from particle physics experiments no longer apply, as the standard 
model fields are confined to the brane. The extra dimensions only have to be smaller 
than the scale on which gravity experiments have probed, of order 1mm at the time of 
writing. Another way to recover four-dimensional gravity at large distances is to embed a 
positive tension 3-brane into an AdSs bulk [9,10]. In this scenario four-dimensional gravity 
is obtained at scales larger than the AdS radius. Randall and Sundrum showed that this 
could produce sensible gravity even if the extra dimension was not compact. 

The cosmology of these extra dimension scenarios has been investigated and the Fried- 
man equation derived and shown to contain important deviations from the usual 4-dimen- 
sional case [11-13]. Some inflationary models have been investigated [14], as have brane 
world phase transitions, topological defects and baryogenesis [15]. 

The possibility that various cosmological problems could be solved in an extra-dimen- 
sional scenario has been examined by several authors; however a novel suggestion that 
could possibly resolve the horizon problem was made in [16], where Chung and Preese 
used a variety of simple metrics to demonstrate that gravity signals, propagating either 
purely in the extra dimension or via a second hidden brane, could connect regions of our 
four- dimensional universe naively thought to be causally unconnected. Such a mechanism 
would have important cosmological consequences as it could greatly alter the size of the 
particle horizon as well as possibly providing a method of experimentally confirming the 
existence of extra dimensions. 

The aim of this present work is to fully investigate both the one and two brane mecha- 
nisms proposed by Chung and Preese, in a much more general, physically acceptable brane 
world scenario: that of a non-Z2 symmetric cosmological model with a non-zero Weyl tensor 
component [17-22]. 

In the single brane case we examine numerically the nature of null geodesies leaving our 
3-brane at various times in the universe: whether they return to our 3-brane or freeze out 
at the horizon; their apparent speed relative to that of light as viewed by a four- dimensional 
observer; and how these attributes depend on the initial velocity of the geodesies. The effect 
on the geodesies of breaking the Z2 symmetry and of varying the Weyl tensor component is 
also investigated fully. We have therefore considerably generalised the previous work done 
on the subject [23-25], showing that there are only significant deviations from Lorentz 
invariance at early times. 

The second case, where signals travel along a hidden brane, has not previously been 
investigated in any realistic scenario: the metrics used in [16] were not solutions to Einstein's 
equations. In order to examine the two brane case fully, the non-linear equation of motion 
of the inter-brane distance (otherwise known as the cosmological radion TZ [26]) is concisely 
derived and is used to evaluate the relative speeds of signals on the second brane compared 
to the first, quantitatively for the stationary case and qualitatively for the moving brane 
case. It is found that although signals can in some situations propagate along the second 
brane significantly faster than along ours, it is not a large enough effect to solve the horizon 
problem; on the other hand we do get significant violation of Lorentz invariance at late 



times such that the gravitational horizon is larger than the particle horizon. This could 
lead to signals in future gravitational wave detector experiments. 

The paper is organised as follows: the general setup and geodesic equations are derived 
in section ^; section |3| contains a brief discussion of geodesies in the normal Randall- 
Sundrum model [23]. We start, in section ^, by investigating non-cosmological (i.e. static) 
solutions of Einstein's equations, with the Randall-Sundrum model extended by admit- 
ting [19] the possibility of Z2-symmetry breaking across the brane, as well as perturbations 
from perfect tuning; the behaviour of null geodesies in the five- dimensional theory is exam- 
ined both analytically and numerically. We then extend the investigation in section ^ to 
fully cosmological solutions of Einstein's equations without mirror symmetry, numerically 
investigating the possibilities for solving the horizon problem in a similar fashion to [25]. 
In section [oj we consider the possible acausality in a two brane scenario where in general 
the second brane is moving. Our conclusions are discussed in section ^. 

2. Geodesies in S-dimensional Einstein Gravity 

Throughout this paper we investigate the behaviour of geodesies in a variety of different 
brane world scenarios. In this section we describe both the general setup and the assump- 
tions that have been made, as well as giving the equations from which all the results in this 
paper are derived. 

In each scenario that we investigate, it is assumed that the five-dimensional spacetime 
satisfies three-dimensional homogeneity and isotropy as we are mainly interested in real- 
istic cosmological models. We choose to work in a 'brane-based' coordinate system as it 
will facilitate physical interpretation of the results. We therefore assume that the four- 
dimensional universe in which we reside is situated on a 3-brane which is chosen to be at 
the origin of the extra dimension (y = 0). This implies that the most general metric must 
have the form 

ds^ = -n2(t, y)dt^ -Fa2(t, y)7y•dxMx^+dy^ (2.1) 

where t is the cosmic time on our brane, x* represents the three spatial dimensions of our 
brane, and y is the coordinate of the extra dimension, jij is the maximally symmetric 
three-dimensional metric with k = —1,0,1 parameterizing the spatial curvature, although 
throughout most of this paper k will be set to zero. The bulk is assumed to be empty 
except for a cosmological constant A, and therefore the metric is obtained by solving the 
five-dimensional Einstein's equations^. 

Gab - ^9AB = fi^TAB, (2.2) 

where we define k^ = I/M5 , M5 being the fundamental (reduced) five- dimensional Planck 
Mass. In the single brane scenario, the energy-momentum tensor Tab takes the form. 

Tab = 5{y)diag{po,Po,Po,Po,0), (2.3) 



''We use the standard brane world convention in that lower-case Roman indices (such as i, j) run across 
the normal space dimensions (1 to 3); Greek indices {fj,, v) run across time and normal space (0 to 3), while 
capital Roman indices {A, B) cover all space and time dimensions (0 to 3 and 5). 



and for the two brane scenario, 

Tab = 5{y)dmg{-po,Po,Po,Po,0) + 5iy-nit))dmg{-p2,P2,P2,P2,0). (2.4) 

where we follow the notation used in [26] by defining po and P2 to be the energy density 
of our brane and of the second brane respectively, and by defining the pressures Pq and P2 
similarly. The position of the second brane given by y = 1^(1), is in general time dependent. 
In Sections (|^ and (^) we use solutions of Einstein's equations that do not possess a 
Z2 symmetry (or mirror symmetry) across the brane. For a full discussion of this topic 
see [17-22]. 

2.1 Geodesic Equations 

Here we derive the necessary geodesic equations corresponding to the metric given by equa- 
tion ( |2.1[) . In what follows the notation is such that dots indicate differentiation with 
respect to the affine parameter and dashes with respect to the fifth dimension, y, leav- 
ing g^ as the differentiation operator with respect to coordinate time. Starting from the 
variational principle 

6S = 6 f£ds = 0, (2.5) 

with 



C = ^OABi^x^, (2.6) 

we can derive the equations of motion for test particles. We shall principally be concerned 
here with lightlike geodesies, which means that we can consider the variation of JC^ds, so 
that the effective Lagrangian density is —n^{t,y)P+a'^{t,y)xiXi+y'^, giving Euler-Lagrange 
equations: 

1 /dn -r, . da 9j9j 



n \dt dt a?n 

Xi = % (2.7) 



o 

y = — '— I ) — y 



2 



where 9i are integration constants. For null geodesies, with which we shall be exclusively 
concerned in this study, we have the additional constraint that the test particle must move 
at the speed of light; this tells us that the first integral of the t equation above must be of 
the form 

nH' = y^+^J^. (2.8) 

These equations will be used extensively throughout sections |3|, § and ^ to evaluate possible 
shortcuts through the bulk. 

3. Geodesies in the Randall-Sundrum Model 

3.1 The Randall-Sundrum Metric 

The Randall-Sundrum model, as initially presented [9,10] is not a cosmological one; rather, 
it is included here as a simple model to develop the reasoning employed in rather more 



complicated cases. It may be derived from the above, general, metric Ansatz given in 
equation (|2.lD by taking n(/, y) = a(/,y) = a{y), and then solving Einstein's equations 
while assuming a Z2 symmetry across each brane. Provided that the bulk cosmological 
constant A is less than zero, and that energy densities of each of the branes are tuned such 
that 

K-po = -K-p2 

the familiar Randall-Sundrum metric is recovered: 



.2„ _ ,.2„ _ ,/ZqX, (3.1) 



ds^ = e-^^ls'lr/^^dx^dx'^ + dy2, (3.2) 

which implies that 

n(y) = a(y) =6-^1^1. (3.3) 



/x is the inverse of the AdSs curvature radius and is given by /i = \/— A/6. 

3.2 RS Geodesies 

We can immediately see that the Randall-Sundrum model is not going to help us solve the 
horizon problem. The local speed of light is everywhere the same, because n{y) = a{y), 
and so even if off-brane null geodesies can return to the brane, they will return with an 
effective sub-luminal velocity. 

However, the Randall-Sundrum model does not exhibit even this kind of behaviour; as 
shown in [23], the geodesic follows the path given by 

g2Ms/ ^ g2/.yo + 2^lyoe^^'yH - v^fih^, (3.4) 

with v^ = 'q^^x'^x'^ . 

As a side note, v^ = does not necessarily correspond to a lightlike geodesic, as claimed 
in section 2 of [23], as the condition for lightlike behaviour is gABX^x^ = 0, which includes 
an extra y^. 

4. Geodesies in Other Static Braneworld Models 

4.1 Obtaining the Metric 

We relax some of the assumptions made in the Randall-Sundrum model; specifically, we 
no longer impose strict Z2 symmetry, nor do we fine-tune the bulk cosmological constant. 
Further to the relaxation of fine-tuning and Z2 symmetry, we must also adapt the Ansatz 
so that n{y) is no longer equal to a{y) to obtain a self-consistent solution. 
The Goo equation from [11] yields 

a^(y) = cosh(2/iy) + Asmh{2fiy), (4.1) 

where we have implicitly scaled so that uq = a\y=o = 1. To consider the loss of Z2 or mirror 
symmetry, when we apply the Israel junction condition 

[a'{y)] = -yP, (4.2) 



where [/] = f{0'^) — /(O ), we make the assumption that a'(O^) = — a'(0 ) + da, giving 

M^ = -yP+y; (4.3) 

da = for Z2-symmetric braneworlds. Then, parameterizing the lack of fine-tuning by 

r = ^ (4.4) 

and the Z2 symmetry breaking by 

we obtain 



a(y) = ^cosh(2^|y|) + (d - r) sinh(2^|y|). (4.6) 

The G55 equation then yields 

^ sinh(2^|y|) + {d-r) cosh(2//|7/|) 

^^ (d - r)Vcosh(2/i|y|) + [d - r) sinh(2/u|y|) ' ^ ' ' 

which has again been scaled so that n{y = 0) = 1; this solution may then be substituted 
into the other Einstein equations and junction conditions to check that it is self-consistent. 

4.2 Geodesies in Static Non-Z2 Brane Worlds 

Since the Lagrangian density £ given by equation ( ^ ) does not depend explicitly on 
coordinate time t, we may obtain a further first integral of the form 

i = -^ (4.8) 

and so straight from the Lagrangian constraint for lightlike geodesies we have 

^0 Gi0i^ 

n^{y) a{y) 



■2 _ ^0 



y = Z^TTT - 777^- (4-9) 



This is effectively a conservation of energy equation, with a kinetic term proportional to y 
and a potential of the form 

ni/) = -^ + 4%- (4-10) 

We may therefore reason that firstly if a geodesic starts its life on our brane (at y = 0) then 
we have, since y^ > 0, 

el > eA. (4.11) 

Further, if the geodesic is ever to return, then y must at some point be zero, so 

„2 \ 



02 M^ Y_ 

°U2(y) a^{y))' 



(4.12) 



where < 7^ = -^ < 1, must be zero at some positive value of y. Rearranging and 
substituting the above solution for the metric, we find that the positions at which y = 0, 
denoted by y^ , satisfies 



or 



+ 1 1-1 



(d-r)(±7-l) 



(4.14) 



{d — r)2 =F 7 

However, there could be a horizon intervening between y = and this point; a horizon 
will occur if either of a{y) or n(y) is equal to zero. We split the analysis into two parts, 
one for (d — r) > —1 and one for (d — r) < —1, reminding the reader that {d — r) = —1 
corresponds to the Randall-Sundrum braneworld. Below, subscript t refers to the geodesic 
turning point, while subscript h to the position of the horizon. 

4.2.1 The {d-r)> -1 Case 

For {d — r) > —1, a{y) (equation |4^ ) is never zero because cosh(2/iy) > sinh(2//y) for all 
y. n{y) = gives (from equation |4.7D 

y>= tanh-^(r-d). (4.15) 



Then y^ (equation 4.14| ) is clearly greater than y^ (as (d— r)^ < 1); also yf is either greater 



than y^ or is negative (the latter when {d — r)^ < 7). Thus a geodesic leaving our brane 
in the positive y direction for (d — r) > —1 will always reach a horizon at y = yj^, and will 
for an observer on the brane therefore never return. See figures [l| and g for illustrations of 
a sample geodesic in this regime. 

4.2.2 The {d - r) < -1 Case 

For (d — r) < —1, n(y) is never zero, and a{y) = gives 

This is always greater than y^~ which is the smaller of the two turning points, and so for 
{d — r) < —1, the geodesic will turn around at a y value of 



yt = TT tanh ^ 



(r-(i)(l-7) 



(4.17) 



(d — r)2 — 7 
Figures ^ and ^ similarly show a sample geodesic for {d — r) < —1. 

5. Cosmological Models 

5.1 The \\0w-Z2 Symmetric Cosmological Metric 

In this Section we investigate the behaviour of geodesies in the full cosmologically realistic, 
non-Z2 symmetric brane world scenario (see [19] for more details). In order to obtain 
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Figure 1: Graph of y against coordinate time for a geodesic in a static untuned {d — r = —0.9) 
braneworld with io = —0.5. The dotted hne represents the calculated y^ horizon position. 



standard cosmology at late times we make the usual assumption, first noted by [12], that 
our brane possesses an energy density which is the sum of the brane tension a and a physical 
energy density p, and we assume that the brane tension and bulk cosmological constant 
are tuned to ensure that the effective four-dimensional cosmological constant is zero. The 
Friedman equation is then given by [19], 



m 



aoV 



K^O 



"1^^+36^ 



k 



+ 



C 



+ 



F' 



(p + a) 



(5.1) 



where ao(t) = a(t, 0), C is the dark radiation term corresponding to a non-zero Weyl tensor 
component, and F represents the extent to which the Z2 symmetry is broken. Motivated 
by previous work on geodesies [24] which suggests that more interesting effects are to be 
seen at earlier times, we assume a radiation dominated universe and therefore set 



A 



ao(t)^ 



(5.2) 



where A is some constant. We can now use the bulk solutions to the metric Ansatz (|2.lD 
which have been found for a single brane without reflection symmetry in an infinite fifth 
dimension [19], and are given by 



«'(i,y) 



Oo 



(i) {A(t) cosh 2/iy + B{t) sinh 2\iy + C(t)) , 



(5.3) 
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Figure 2: Graph of x against coordinate time for a geodesic in a static untuned (d — r = 0.9) 
braneworld with xq — —0.5. Note that as the geodesic goes towards the horizon its x- velocity 
decreases. 



where ^ as before, is defined in terms of the bulk cosmological constant A as ;U = ^— A/6. 
The purely time dependent constants A{t), B{t), and C{t) are given by 



A{t) 
B{t) 
C{t) 



1 + (1 + ti;)- + -^ + 



2^2 



PP 



<y (a + p) 



a 2cj2 2(a + /3)2 
fp 



2„2 



-(! + «;: 



fp 



a 2 a 2((7 + p)2 



(5.4) 
(5.5) 
(5.6) 



Here w and / are dimensionless constants defined by w = 18C/ko"A and / = QF/ tt^aX, and 
the ib-signs in the expression for B{t) give the two different solutions on either side of the 
brane. The solution for n{t, y) is given in terms of a{t, y) as: 



n{t,y) 



do(t) ' 



(5.7) 



where from now the dots denote differentiation with respect to coordinate time. 

To find ao(t), we must solve the Friedman equation ipl\); however, for / / 0, there 
is no analytical solution, and so it must be solved numerically; this presents a problem, 
however, due to the initial singularity. Therefore, to solve this equation we approximate 
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Figure 3: Graph of y against coordinate time for a geodesic in a static untuned (d — r = —1.1) 
braneworld with xq = —0.5. The dotted Hne represents the calculated y turnaround position. 



the /^ term for early times, when p ^> o", as f'^o''^] the solution is then, as in [19], 



ao 



7/ 



fcT 



2af 
cosh — r^-t — 1 
3M| , 



1 . , 2af ' 

-\ smh — ^-t 

cr 3M| ^ 



(5.^ 



This solution will be approximately valid to some time Iq; as an estimate of that time, we 
find when p = Aa, giving 



or 



to 



l + w 
3Mi 



cosh — ^^tn — 1 
. 3M| 



+ smh — ^-to = -J, 
3M| 4 



2af 



log 



^ + { + V(^ + i)'-(^)' + i 



1+"' + 1 



r 



(5.9) 



(5.10) 



1 3Af| 



see 



Over the range of values of interest [0 < w < 0.2, < / < 60), to > 20^2^ 
figure ^) and so we may safely take the initial conditions for numerical integration of the 
Friedman equation as our approximate solution at t = 40^^- This can then be numerically 
integrated using the full Friedman equation ( ^ ) to generate the solutions of the Einstein 
equations. 

5.2 Geodesies 

We would expect at late times for geodesies to be similar to the static case, assuming that 
the dynamical timescale is much smaller than the cosmological one. This appears broadly 
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Figure 4: Graph of x against coordinate time for a geodesic in a static untuned {d — r = —1.1) 
braneworld with xq = —0.5. The dotted hne represents a Hght ray moving on the brane. Note that 
the true nuh geodesic has covered a greater distance on return to the brane than the Ught ray in 
the same time t; lower magnitudes of io lead to higher superluminosities. 



to be the case, though the position or cross-section of the horizon would need to be tuned 
for an adequate comparison. The main difference is that for large dynamical times there 
will be a certain amount of asymmetry in the geodesic's path, due to the changing scale 
factor a = a{t,y). 

Consider figure |6|, which shows null geodesies leaving the brane at a time Iq such that 
ao(^o) = 0-5, when the p'^ term dominates in the Friedman equation^. In that figure, both 
geodesies have an initial y of 0.975 (where the speed of light is 1); the difference is that 
the dotted geodesic is in a braneworld with broken Z2 symmetry (/ = 0.1, in the notation 
of [19]). The graph has been cut off so that the path of the returning geodesic is clear; the 
/ = 0.1 geodesic asymptotically (as t tends to infinity) reaches the horizon in the bulk. It is 
perhaps slightly clearer when viewed in conjunction with figure 0, which shows that when 
the geodesic in the / = 0.1 world reaches x of about 290 it is at the horizon, at which point 
it will move with the horizon in the y direction. 

Let us next examine figures § and |9|, in which we have taken a number of geodesies as in 
figures ^ and 0, and examined the superluminosity of the returning signal; in other words, 
how much faster than light, which is constrained to travel along the brane with speed ^^^jy, 
is the true null geodesic. 

The results are perhaps surprising at first sight. Conceptually there are two separate 



^Though extensive tests have shown that there is no particular difference in between the regimes. 
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Figure 5: A graph of to (for w = 0) over a range of values of /; this is used to determine the point 
at which numerical integration of the Friedman equation must begin. 



effects happening here. One is that space is 'warped'; loosely, there is an e'^^ factor warping 
space in the bulk, meaning that null geodesies can cover effectively much greater distances. 
Thus, the farther into the bulk that the geodesic penetrates, the more superluminal the 
geodesic will be. However, there is a competing effect due to the presence of a horizon in 
the bulk; the horizon is only a coordinate singularity, but since it is a singularity in the 
physical brane coordinates it is nevertheless of physical relevance, in that geodesies that 
reach the horizon never return from the point of view of an observer on the brane. 

Null geodesies that have too large a velocity component in the fifth dimension will not 
be able to escape the horizon. Thus, at very high yo, the geodesies will not return to the 
brane before the end of the universe. Null geodesies with small initial y, however, will be 
confined to the brane. 

In between these two regimes, what will happen? At some critical angle, the effect of 
the horizon will begin to dominate over the effect of the warped spacetime, and consequently 
there will be a maximum in the ratio of gravitational to electromagnetic effective speeds. 

As discussed in [24], the maximal superluminal effect will be obtained when the hier- 
archy is the most pronounced experiment ally- allowed value, and the null geodesic is set off 
at the earliest physical time possible (of the order of M^ ). However, it is not because of 
any particularly stronger warping at this earlier time, but rather because the geodesic can 
spend much 'longer' in the warped area before returning to the brane. Thus there is no 
discrepancy in principle^ between the relatively low superluminal values shown here and 



^Unfortunately it is not possible without extensive code writing to numerically investigate this figure in 
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Figure 6: Graph of y against coordinate time for a geodesic in a cosmological braneworld with 
yo = 0.975. The returning geodesic is in a Z2-symmetric braneworld, while the one tending towards 
the horizon has / = 0.1. The geodesies were set off at oq = 0.5 



the figure of 10^ quoted in [24]. 

It appears from figures ^ and ^ that the effect of including asymmetry^ is to increase 
the horizon cross-section, for geodesies leaving at equivalent times (taken as equal redshifts, 
or, what amounts to the same thing, equal energy densities), and also to decrease the ratio 
of the speeds even for geodesies that do not go anywhere near the horizon. 

The effect of the start time on the superluminosity is such that when the p^ term 
dominates, at early times the superluminosity is increased (see figure |^). A corresponding 
late-time graph is not shown, as the cross-section of the horizon has increased to such an 
extent that almost all geodesies hit the horizon and do not return, from which we note that 
the acausality in this model is maximal at early time, with energy leakage from the brane 
at later times. 



Including dark radiation likewise dampens the superluminosity; figure |11| shows the 
effect of including a dark radiation term tt; = 2 in the calculation. Note that this value of 
w is well outside the nucleosynthesis bound [27], and is tested in our investigation only so 
as to exaggerate any possible visible effect. 



this framework, as the maximal superluminosity will be obtained with an initial ij of somewhere between 
1 - lO"^'^ and 1. 

'^We have shown here results with positive / only, as results with / < are physically similar. 
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Figure 7: Graph of x against coordinate time for a geodesic in a cosmological braneworld with 
yo = 0.975. The geodesies are the same as in Figure || Note that the / = 0.1 geodesic asymptotically 
tends to a maximal x = x{th), where th is the value of the afHne parameter when the geodesic 
crosses the horizon. 



6. The Two Brane Model 

It would seem from the above results, that the possibility of null signals taking shortcuts 
through the bulk in the single brane scenario would not solve the well known cosmological 
horizon problem as the apparent speeds of such gravity signals are not significantly greater 
than the speed of light signals confined to the brane. Due to this, we now investigate an 
alternative suggestion made by [16], which was that a 2-brane scenario could provide a 
solution to the horizon problem. By investigating a variety of 'toy' metrics, Chung and 
Freese showed that apparently acausal signals could be sent between different points on 
our brane, via the second hidden brane. These signals would, as before, bring in to causal 
contact regions of our Universe that would have otherwise been unable to communicate. 
In order for this method to work, a metric of the following form was assumed: 



ds^ 



dt" 



-^^'ya^{t)dx^ 



dy^ 



(6.1) 



where y corresponds to one extra spatial dimension. Our brane and the hidden brane were 
assumed to be at y = and y = TZ respectively. The distance Dae travelled by a null 
signal on our brane in time tj is then (see figure p^ ) 



Dae 



'S dt 



a{t) 



(6.2) 
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Figure 8: The ratio of gravitational to electromagnetic speeds, as a function of the initial compo- 
nent of velocity in the fifth dimension, in a Z2-symmetric braneworld. ao(io) = 0-5. 



However, if it is possible for a null signal to leave our brane, interact with and therefore 
travel along the hidden brane, before eventually returning to our world, we can then ask 
how far it would have appeared to have travelled Dad, in time tf. If it takes a time tc to 
cross between the branes then this distance would be given by 



Dad 



oi^T^ 



a{ty 



(6.3) 



Therefore if e^ ^ 1 then Dad 3> Dae and the horizon problem will be solved. The major 
problem with this suggestion is that it involves metrics that are not realistic. It relies upon 
there being a conformal factor that affects the spatial part of the four-dimensional metric 
only. In typical Randall-Sundrum Brane World models the metric is usually of the form 



ds^ = e-^^yri^^dx^dx" -dy^, 



(6.4) 



and hence the horizon problem will remain unsolved. 



6.1 A Cosmologically Realistic Metric 

are obviously not realistic descriptions of our Universe; 



The above metrics (|6.l| ) and 

it is therefore worthwhile to consider if the above proposed solution to the horizon problem 

will work for a cosmological metric. liw, f and k are set to zero, the solutions to Einstein's 

equations given by ( ^.3| ) can be greatly simplified. Writing a{t, y) in terms of exponentials 

gives 



a{t,yf 



ale-^f^y 



l^2^y _ (£ + 2) 
a a 



(6.5) 
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Figure 9: The ratio of gravitational to electromagnetic speeds, as a function of the initial compo- 
nent of velocity in the fifth dimension, in a non-Z2-symmetric braneworld (/ = 0.1). Note firstly 
that the maximum has moved leftwards, and also that some geodesies from our samples no longer 
return. 



This can be further simplified to give 

a{t,y) = ao (cosh /iy - r/o sinh ;uy) 
which can be used with equation ( |5.7D to give 

n{t,y) = {cosh ij,y - fjo sinh fiy) , 
We have defined 



Vo 



1 + ^, 
a 



m 



^o + jr- 



(6.6) 



(6.7) 



(6.5 



Here, Hq is the Hubble constant on the brane at y = 0, given by equation 5.1 and can now 
be written in the form 



Hi 



K^ 2 ^ 

36^° "6 



^ 



(*= - 1) . 



(6.9) 



In order for there to be a significant difference in the relative 'speeds' of null signals travelling 
along the two branes, we require that n/a is larger for some y than it is for y = 0, i.e. that 



n{t,y) n{t,0) 

—, r 3> — — -7 -^ 

a{t,y) a{t,0) 



cosh fiy — fjQ sinh ^y 1 

ao(cosh ^y — rjQ sinh /iy) oq 



(6.10) 



An acceptable equation of state on our brane demands that p oc Oq "^ and this implies that 
fjo = Vo ~ QP/^ ^-iid therefore that fjo < rjQ. Because of this the ratio n/a will diverge as 
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Figure 10: The ratio of gravitational to electromagnetic speeds, as a function of the initial com- 
ponent of velocity in the fifth dimension, in a ^2-symmetric braneworld. ao(^o) = 0.125 (earlier 
time than Figure 0) . 



tanh iiyh 



and is given by 

(6.11) 



y -^ yhi where yn is the position of the horizon defined by a{t,yh) 

1 

m' 

This suggests that null signals may travel along the second brane at much greater 'speeds' 
than on our brane, provided the second brane is close to the horizon. Unfortunately, to 
correctly calculate the difference in speeds, we need to know the motion of the second brane 
with respect to the one at y = 0. 

6.2 General Behaviour of the Inter-brane Distance 

In this section we use a simple non-perturbative method to derive the general equation of 
motion for the inter-brane distance otherwise known as the cosmological radion and denoted 
by 'R.{t), which was first derived by [26]. For alternative approaches to this topic see [28-31]. 
Noting first that if the second brane follows a trajectory given by y = 7?.(t), then the 
induced 4-dimensional metric on the second brane is given by 



ds' 



n^{t,n{t))-V? dt^ + a^{t,n{t))dx^ 



-dr"^ + al{T)dx^ , 



(6.12) 
(6.13) 



where r has been defined as the proper time as seen by an observer on the second brane. 
The expansion rate of the second brane as seen by an observer on our brane is simply 

02 di \a a / 2 



H2{t) 



(6.14) 
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Figure 11: The ratio of gravitational to electromagnetic speeds, as a function of the initial com- 
ponent of velocity in the fifth dimension, in a Z2-symmetric braneworld. w = 2.0 (same time as 
Figure 0). Note that some previously-returning geodesies no longer return. 
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Figure 12: A Null signal that remains on our brane goes from A to E in time tf, one which leaves 
the brane could travel between A, B, C and then return to our brane at D in the same time as the 
first signal reaches E. 



and as seen by an observer on the second brane itself 

1 da2 



W2(r) 



02 dr 



- + -n)( 

a a J ^ 



-<r 



^n^ -V} 



-1/2 



(6.15) 
(6.16) 



This equation for 7^2 is important as it relates the expansion rate of the second brane (with 
respect to proper time) to 7^ and 7^. It does not seem so useful at first, as calculating 7^2 
could be difficult; however, we know that the brane world Friedman equation given by ( |6.9D 
is derived from a purely local analysis and that should we have chosen the second brane to 
be stationary and at y = 0, we would have derived the equivalent Friedman equation with 
p replaced by p2- This means that 7^2 ("T") must have the form 

nlir) = ^^pI-^ = ^?{r,l-l). (6.17) 

Where we have defined r]2 = P2/o'. Equation (|6.17|) ensures that the second brane is Z2 
symmetric and that it evolves according to the junction conditions of the extrinsic curvature 
tensor, just as equation ( |6^ ) ensures similar behaviour for our brane. At this point we use 
the following identity^ which is obtained from the Einstein equations and is derived in [27]: 



a \' a'2 



-V 



naj a? 



^i\ (6.18) 



givmg 



m = (£-)'- ^ + Ai. (6.19) 



naJ a^ 



Substituting ( |6.19 ) into ( |6.15| ) and rearranging finally gives the following first order equation 
for 7t: 



f^. , \ / a2 



-1 



^ = ^("S ^^'^^V^^)ul^ + ^''?2) ■ (6.20) 



This can be rewritten in the form found by [26] 



a' a ■ 




( 


7?,^ 


1/2 


- + -^n = 


= f^V2 


i- 






a n'^a 




\ 


'/r 





(6.21) 



In order to investigate acausal signals we will be interested in a stationary second brane, and 
therefore need an equation for TZ. [26] have shown that by differentiating ( 6.21] ) with respect 



to time, considering the three-dimensional symmetries of the bulk energy momentum tensor, 
and the Z2-symmetry of the branes, it is possible to derive the following equation for TZ 



Equations ( |6.20 ) and (|6.22D govern the evolution of IZ and will be used in the next two 



sections to evaluate the difference in the speeds of null signals travelling along each brane. 



Note that although the derivation is done here for / = fc = ui = 0, this method generalizes to 
braneworlds with non-zero values of these para meters with the use of the appropriate Friedman equation 
for equation |6.17| and identity for equation |6.1q. 
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6.3 Acausal Signals: Stationary Branes 

The simplest solutions to equations ( |6.20| ) and ( 6.22 ) are those describing a stationary 



second brane; in this section we will investigate the appropriate conditions and the ratio of 



distances travelled on each brane in a certain time interval. Setting 7^ = 7^ = in (|6.2lD 



and ( 6.22| ) trivially gives the well known conditions on r/2 and p2 = i-^/c which ensure that 



the second brane does not move: 

^ sinh /x7^ - 7?o cosh /i7^ 
cosh iJiR, — rjo sinh fiTZ ' 
1 sinh fiTZ — fJQ cosh fiTZ 2 
3 cosh fiTZ — fjQ sinh ^TZ 3 

Note that if we demand standard cosmology on our brane at late times then this requires 
that both r/o and f/o tend to 1 (not zero as [26] have suggested). This results in the 
equation of state of the second brane becoming — r/2 = P2 = 1 which is equivalent to 
—p2 = P2 = 6/i/K^, the conditions in the original Randall-Sundrum two brane scenario. 
Assuming the above form for the energy and pressure densities on the second brane, we 
now just have to evaluate the distances travelled by null signals between times ti and ^2 on 
either brane: 



Di = / —At, D2 = I ' ' ;,;:,, dt, (6.25) 



hi ao ' Jti a?{t,TZ{t)) 

however since we have assumed a stationary second brane, D2 simplifies to become 



D2 = f!^MMdt = r cosh^7^-.)osinh^7^ ^^^ 
Jti a{t,TZ{t)) Jti ao{cosh fiTZ — rjo sinh fiJl) 

As mentioned before, a standard equation of state on our brane leads to fjo having the form 
% = % — QP/(^ and therefore. 



D2 = / 

Jti 



1 {qp/a) sinh pR, 



ao ao(cosh//7^ — (1 + p/o") sinh/i7^)J 
= Di + AD. (6.28) 

Rewriting the Friedman equation for our brane in terms of p/a gives Hq = p^ q^ {{p / (y)"^ + 
2p/a). This, and the fact that p = A/oq allows the IS.D integral to be converted into the 
following form: 

piaXfli Jm ui/<?^/2^rfT(n - A{n)) ^^''^^' 

where u = a/p and all the information on the second brane's position has been grouped into 



A(TZ) = sinh ^7^/ (cosh ^7^ — sinh//7?.). The integral in ( 6.25 ) can be similarly converted, 
and we then obtain the ratio of the extra distance travelled on the second brane, divided 
by the distance travelled on our brane in the same time, 

AD ( /""^ dti \ f /""^ d-u \ ^^ 

— = qA{n) U^ -^-^==-——— ^^ ____j . (6.30) 
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We now need to determine the possible limits on ui and U2- Obviously the contribution 
of the integral to AD for values of u greater than {q + 1)74(7^) will be less than the cor- 
responding contribution to Di and so we choose the upper limit on U2 to be {q + 1)^(7^), 
as a larger U2 is of no interest. The minimum possible value of Ui is slightly harder to find 
as it is constrained by requiring that the energy densities of both branes remain below the 
5-dimensional Planck mass limit 

_6_ 

M' 



< 



Mi 



a 

- < 

P 



6M| 



(6.31) 



where we have defined M to be the dimensionless ratio of the squares of the four- and 
five-dimensional Planck masses M = M^/M^. Demanding a similar constraint on p2 and 
using the stability condition given by ( |6.23|) , results in the following condition 

I no I M 

(6.32) 



\P2\ 



< 



Mi 



\V2\ 



a 
r/o cosh nTZ 



< — 

- 6 

sinh iiTZ 



cosh fj,TZ — rjQ sinh fj,TZ 6 

This can be rearranged and expressed in terms of u using the fact that r]Q = 1 + 1/u, to 
give the constraint on u that ensures that \p2\ is always less than M5 : 



> 



M tanh fill + 6 



(M-6)(l-tanh^7^)■ ^^'^^^ 

In the next two sections we will use this constraint to determine ui which corresponds to 
the earliest possible time that a null signal can set off along the second brane. We examine 
both the near brane (pTZ <^ 1) and the far brane (pTZ ^ 1) limits. 



6.4 Near Brane Limit 



We will now evaluate the integral expression for AD/Di given by ( 6.30 ), in the small fiTZ 
limit. The minimum possible value of the inter-brane distance 7^ will be the inverse of 
the 5D Planck mass: 7?. > I/M5 and for a cosmologically realistic brane fi is given by 
// = M^/M^, leading to the relation 



fxTZ 



> 



Ml 

Mi 



1 

M' 



(6.35) 



If we therefore assume 1/M < fiTZ ^ 1, we then find that in this limit A{TZ) ~ fiTZ and 
therefore that ui < U2 = {q + 1)^(7^) ^ 1. Equation ( 6.30| ) can then be approximated by 



AD 
1^7 



qA{n) / 

\Jui 



{q+l)A{n) 



du 




(g+l)A(7^) ^^ 



-1 



aV? 



If we now assume that our brane is undergoing radiation dominance by setting q 
first integral on the right hand side of ( |6.36 ) can be integrated to give 

du 



(6.36) 



4, the 



Al 



Ul 



u^/'^iu - A{n)) 



2 arctan 



2 arctan 



+ In 



(^)3_l (^)i+l 



(^)3_1 (^)3+l^ 
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and the second integral is trivially given by 



c\i/ 433 

-{u-^-ut). (6.37) 



m nV4 3 

From above we take U2 = {q + 1)A{TZ) = 5A{TZ), and therefore we just need to evaluate 
ui/A{TZ) in the small /iTZ limit using equation ( |6.34D 

ui Mtanhu7^ + 6 3 , , 

A{R) (M - 6)tanh/x7^ 2M^7^■ ^ ' ' 

Replacing the values for ui and M2 into the expression for IS.D /Di and ignoring all sub- 
dominant terms, finally gives the ratio of the maximum extra distance travelled by a null 
signal on the second brane compared to the first in the near brane limit as 

— - ~ 1.281n(M^7^). (6.39) 

The value of M5 is constrained by nucleosynthesis such that M5 > 30TeV. This leads to 
the corresponding constraint on the dimensionless ratio M to be 

M = ^ < 10^9. (6.40) 

In order for the above approximations to be valid ^TZ must satisfy iilZ < 10"^, therefore 



(6.39) becomes 

AD 

— - ~ 80. (6.41) 

This shows that any seemingly acausal effects due to signals travelling on the second brane 
as opposed to our brane are too small to solve the cosmological horizon problem in the 
(7^ <C 1/^) limit. Another important point is that /S.D is only significantly greater than 
Di for a very brief time, shortly after the big bang. If the signals were to travel for 
longer periods, say until ti2 = 1 then the distance ratio would become much less than 1: 



AD/Di ~ iJ.TZlia{M fiTZ) as can be seen from ( 6.36 ). 

6.5 Far Brane Limit 

In this Section we investigate the behaviour of null signals travelling on the second brane 
when the inter-brane distance is large (7^ » 1/m)- The fact that the second brane has to 
be closer to our brane than the horizon implies that we will be examining signals that are 
travelling at late times with respect to our Universe, as opposed to the previous Section 
where the only interesting situations occurred at early times. 

In the TZ ^ 1/fi limit, the function A{TZ) becomes very large and can be approximated 
by 

^ ^ l-tanh/i7^ 2 ^ ' 

As before, we will require the ratios U2/A{TZ) = {q + 1) and ui/A{TZ) which is obtained 



from equation (|6.34j) and is approximately given by 



ni Mtanhiin + 6 ,12 , , , 

"■ 1 + 77- (6-43) 



A{R) (M - 6)tanh/i7^ M 
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Therefore both ui and U2 are exponentially large and so the expression for AD/Di given 
by equation (|6.30|) can be simplified in the large 7^ limit to 



AD / /"W+i. 



iq+l]Ain) 



du 




iq+l)A(n) 



(l+fM(7^) 




(6.44) 



i{i+j^)Ain) u-^ + 2 (^u - A{n)) ^ 

Assuming radiation dominance, the first integral on the right hand side of ( 6.44| ) can be 
solved to give 



Ai 



U2 



du 



U\ 



u3/4(u-A(7^)) 



-2arctan 



+ 2 arctan 



+ In 



(1)^ 


-1 


(1)^ 


-1 



(1)^ + 1 
(1)^ + 1. 



which when combined with ( 6.44| ) and ( |6.43 ), and after all subdominant terms are neglected, 

results in the following expression for the extra distance travelled on the second brane in 

the far brane limit 

-— ~ 21n(M). (6.45) 

Di 

Note that if instead of radiation domination we had chosen matter domination and set 

g = 3, the leading term given in equation ( 6.45| ) would not have been altered. Replacing 

into ( |6.45 ) the maximum possible value of M ~ 10^^ leads to the final result 

— - ~ 130. (6.46) 



Unlike when IZ <C 1//U, here AD is significantly greater than Di for an extended period 
of time. Solving the Friedman equation gives the late time relation between u and the cosmic 
time t as measured on our brane to be u ~ /i^Q'^t^/2. This implies that AD > Di during a 
time interval of At ~ (M4/2MI) exp(/i7^) = 10^'^exp(//7^)TeV~"^. However, as can be seen 



from equation 5.44 , if limits for the integral for AD jDx are such that M2 3> iti, then any 
acausal effect becomes negligible. Taking u\ to be just before recombination and U2 at the 
present day leads to u\ ~ 10^^ and U2 ~ 10^^ , and therefore the horizon problem cannot be 
solved in this manner, though this effect may have important consequences in other areas 
of cosmology such as structure formation. One should note that the effectiveness of the 
acausality would also be lessened by the time taken for signals to cross between the branes. 

6.6 Acausal Signals: Moving Branes 

Up till now we have been examining the possibility of acausal signals travelling along sta- 
tionary branes. It is however interesting to consider the effect of allowing the branes to 
move. As will be shown, for a general cosmological scenario the equation of motion of 
the second brane is very difficult to solve; various qualitative features can, however, be 
investigated and will be discussed in this section. We start with the equation for the time 
dependent inter-brane spacing 7^(t), derived previously ( |6.21 ) 



n 

n 



a' a 



a?n 



+ f^ V2 



V2 



2„2 + /^'^2 



n^a 



(6.47) 
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This can be rewritten using equations (|6.17| ), (|6.ig| ), ( |6^ ) and ( |6^ ) as 

7t HQ{i]QCOshfiTl — sm]ifj,TZ)+firi2TC2{coshfiTZ — r]osinhfiTZ)'^ 

n {rjQ cosh fiTZ — sinh fiTZ)"^ + Ti2{cosh. fiTZ — tjq sinh fiTZy 

For the trivial case where 772 is constant and therefore 7i2 = 0, equation ( |6.48| ) can easily 
be solved numerically as was done in [26]. For a cosmologically realistic case where r]2 is 
time dependent, the solution is less easily obtained. If, for example, we assume a standard 
equation of state on the second brane: P2 = UJ2P2 which provided 0)2 > —1, implies that 
P2 oc l/a2^, we can then solve the Friedman equation 



2 _ i aa2 



da2\' 

dr ) qlril 




T<i = -A^] = -^2\^\ =f^'ivi-l), (6.49) 



to give r]2 = —l/sm{q2HT) and TC2 = /i^/tan(g2/ur). However, to solve equation ( 6.48 ) we 
need to know r/2 and 'H2 in terms of t: our cosmic time, which is related to the second 
branes cosmic time r by 



dr = yn'^-TZ'^dt. (6.50) 

This makes it much more difficult to solve for TZ and due to this added level of complication 
we leave a more general study of ( |6.48D to future work. We now examine the situation 
qualitatively and argue that acausality in the moving brane model is not appreciably greater 
than in the static case. 

From equation ( |6.47 ) it can be seen immediately that since a' < for < i? < y/i, then 



if 772 > then 7^ > for all R and the second brane will move away from our brane and 
freeze out at the horizon, which has a time dependent position given by yh = arctan(l/r7o). 
If ?72 < 0, then 7^ = only when 

Requiring also that TZ = gives the conditions for stationary branes as discussed 
previously. The stability of a stationary second brane has been investigated in [26] where 
it was shown that if both branes are de Sitter the equilibrium position is unstable, while if 
both branes are anti-de Sitter then it is stable. For our purposes, we only need to examine 
the case where the second brane moves in the positive y direction and freezes out at the 
horizon, as the other possibilities are that the second brane approaches our brane and 
collides with it which is physically unacceptable, or that the second brane is stable: the 
case which has been examined previously. 

If the second brane does freeze out, then TZ ^ yh and 7?. — > n{t, R) as can be checked 



from equation ( [6.47[ ). The expression for the distance travelled on the second moving brane 

is 



rt2 ^/n2(i,7^(t))-7^2 

which will therefore tend to zero as 7^ — > n{t,R). Another problem is the distance that 
signals have to travel to return. In order for there to be a significant acausal effect, the 
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null signal has to be travelling on the second brane for a substantial amount of time and 
therefore would be a large distance away from our brane when starting to return. The 
combination of these two effects: the distance travelled by a null signal on the second brane 
frozen at the horizon tending to zero, and the large distance returning signals would have 
to travel, suggest that little or no acausal effects would be observed on our brane. 

7. Discussion 

We have investigated the behaviour of null geodesies in several five-dimensional brane world 
scenarios. In the single brane case it was shown that apparent causality violations caused 
by such signals taking shortcuts through the bulk were small. The ratio of the speed of 
five-dimensional gravity signals to the speed of four- dimensional light Cgjc^ was found to be 
in general not more than two. The effects of introducing a non-zero Weyl tensor component 
and of relaxing the mirror (or Z2) symmetry were examined and shown in general to decrease 
the observed acausality. 

It was found that during early times in the Universe many geodesies would return to 
our brane, as opposed to late times when most of the geodesies leaving our brane freeze out 
at the horizon. This is understandable, as at late times {p -^ 0) the cosmological metrics 
that were considered tend to the standard Randall- Sundrum metric which, as previously 
discussed, results in geodesies not returning. The behaviour of these geodesies could be 
used in conjunction with the temperature dependent production rate of gravitons to exactly 
determine the amount of energy lost to the bulk throughout the history of the Universe as 
is discussed in [32]. 

Analysis of the two brane scenario (with w = f = k = 0) shows that signals can in 
some situations travel along the second brane significantly faster than along our brane. 
This effect, however, would only last for a certain time depending on the position of the 
second brane. In the near brane limit (pTZ <C 1) this period is very brief and any acausal 
effects would become negligible when longer time intervals are considered. In the far brane 
limit (//7^ 2> 1) however, the period of interest can last for At ~ 10^^exp(^7?.)TeV~^and 
it was found that the apparent speeds of null signals were approximately 130 times faster 
on the second brane than on the first. It should be noted that a non-zero w and / would 
lessen the above effect as the crucial ratio n{t,TZ)/a{t,TZ) is largest only when w = f = 0. 
Another factor that would detract from the effectiveness of this mechanism in the far brane 
limit is the time taken for signals to travel between the two branes. 

In neither one nor two brane scenarios are the possible acausal effects significant enough 
to enable us to solve the cosmological horizon problem. The increase in the gravitational 
particle horizon could, however, have important effects on other areas of cosmology such as 
the initial conditions of inflation, structure formation and the production of topological de- 
fects. It would also technically be possible to measure the time delay between the detection 
of gravitational waves and light waves from a particular cosmic event; however, we believe 
that this would require significant advances in gravitational wave detector technology. 

There is also the possibility that a higher dimensional model could produce a large 
enough acausal effect to solve the horizon problem [33-40]. Unfortunately, we are currently 
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unaware of any metrics in {d+ 1) dimensions with d > 4 where acausal signals are possible. 
For example there have been several six-dimensional 'brane world' models proposed that 
have a metric Ansatz of the general form 



ds' 



<i)^{x')T]af3{x'')dx"dx'^ + gjk{x')dx^dx'' 



(7.1) 



where a, /3, and v run across time and normal space dimensions (0 to 3) and i and j run 
across the two extra spatial dimensions (5 to 6). It can be seen that no acausal effects are 
possible since both the time and the three-dimensional spatial components of the metric 
have the same dependence on the extra dimensions, unlike the five- dimensional case where 
^(^,1/) / o.{t,y)- A fully cosmological six-dimensional metric could however overcome this 
problem. 
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